Equilibrium Properties of Trapped Dipolar Gases at Finite Temperatures 



Yuki Endo^, Takahiko Miyakawa^ and Tetsuro Nikuni^ 
^Department of Physics, Faculty of Science, Tokyo University of Science, 
1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan 
'^Faculty of Education, Aichi University of Education, 
1 Hirosawa, Igaya-cho, Kariya, Aichi 44-8-8542, Japan 
(Dated: September 14, 2010) 

We study the equilibrium properties of dipolar Bose and Fermi gases at finite temperatures. We 
recently developed a variational ansatz for the phase-space distribution function of a dipolar Fermi 
gas at finite temperatures. We extend the variational formalism to a Bose gas and discuss the effect 
of dipolar interactions with different statistics on the thermal equilibrium, with particular emphasis 
on the deformation in momentum space. We examine the stability of the system by varying the 
temperature, trap aspect ratio and the dipole moment. In addition, we discuss how deformation in 
both real and momentum space can be observed in the high-temperature regime, which is relevant 
for the current experiments at JILA [Ni et al. Science 322, 231 (2008)] and at University of Tokyo 
[Aikawa et al. New. J. Phys. 11, 055035 (2009)]. 

Interest in dipolar gases has been growing since Bose-Einstein condensates (BECs) of ^^Cr atoms, which have large 
magnetic dipole moments, were experimentally observed [l|-|3|. The anisotropic and long-range nature of the dipolar 
interaction confers interesting properties to the equilibria and dynamics of dipolar gases. A number of experiments on 
^^Cr BECs have investigated their ground states and expansion [l|, [H, 0, HI- On the other hand, several experiments 
are proceeding energetically toward creating heteronuclear polar molecules whose large electric dipole moments give 
rise to strong dipolar interactions. An example of a dipolar Fermi molecule is ^^Rb^^K Q and examples of dipolar 
Bose molecules are ^^Rb^^K [7], ^^^Cs^Li [8], and ^Li^^K [9|. Rydberg atoms have also been attracting interest due 
to their large electric dipole moments. 

There have also been several theoretical studies of dipolar BECs, which have investigated their ground states p^ol^fTlj, 
collective oscillations [l2|, IHI , and their properties in optical lattice potentials [13 , fTHj . There is also growing interest 
in dipolar Fermi gases. The ground state [l6|, [l^ and stability fl?!, Iisj of dipolar Fermi gases has been theoretically 
studied. In addition, there have been studies of free expansion [19, 20], collective oscillations and the possibility 
of superfiuid phases (22l-[25|. 

We are currently focusing on experiments with ^^Cr BECs. For example, the expansion dynamics of ^^Cr BECs have 
been observed and it was found that the long-range and anisotropic character of the dipolar interaction leads to an 
anisotropic deformation of expanding ^^Cr BECs [ll,[2[. These observations indicate that the dipolar interaction tends 
to elongate dipolar BECs in the dipole moment direction and shrink it radially. The aspect ratio of dipolar atomic 
gases after free expansion is very sensitive to the trap geometry and the ratio between the dipolar and interatomic 
interactions. The reason why dipolar BECs elongate along the dipole moment, where two individual dipoles attract 
each other, can be understood by considering the shape of the dipole-dipole potential 0. It was found that the 
dipole-dipole potential within a cloud of atoms is a parabolic saddle with a negative curvature along the dipole 
moment. Consequently, the total energy of a trapped dipolar BEC is reduced if the atoms are redistributed from the 
repulsive to the attractive direction, even at the cost of increasing the external trapping potential energy. 

As mentioned above, the shape of dipolar BECs strongly depends on the strength of the effective dipolar interaction. 
Experimentally, one can adjust the strength of effective dipolar interaction by using a Feshbach resonance, which can 
control the interatomic interaction. By reducing the interatomic interaction, one can enhance the effective dipolar 
interaction. The aspect ratio of a dipolar BEC after free expansion has been experimentally observed using this 
technique. [^]. In this experiment, the large deformation induced by the dipolar interaction was observed by adjusting 
the scattering length to approximately zero, where the effective dipolar interaction is largest. On the other hand, 
many theoretical studies have shown that a dipolar BEC becomes more elongated in the dipole moment direction with 
increasing strength of the effective dipolar interaction [lo|, [O, [HI- In terms of the stability of dipolar BECs, it was 
found that the critical scattering length strongly depends on the trap aspect ratio. In a cigar-shaped trap, a repulsive 
interatomic interaction is required to keep the system stable, whereas in an oblate trap, there is stable region with an 
attractive interatomic interaction because of the effective repulsive dipolar interaction [5[. 

On the other hand, there is a growing expectation of creating dipolar Fermi gases. Unlike dipolar BECs, single- 
component dipolar Fermi gases do not interact via 5-wave collisions and they have Hartree direct and Fock exchange 
energies of the dipolar interaction in the mean-field description, which reflects the antisymmetric many-body wave 
function. Miyakawa et al. introduced a variational Wigner function to examine the ground state of the system at zero 
temperature and they demonstrated that the anisotropic nature of the interaction causes Fermi surface deformation 
through the Fock exchange energy [l^. They found that the real-space distribution, which depends on the trap 
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geometry, is elongated more along the dipole moment than the real-space distribution of system without a dipolar 
interaction and that the shift is largest in an isotropic trap. In addition, they found that the momentum distribution 
is always elongated in the dipole moment direction by the Fock exchange energy. Their variational calculations are 
highly suitable for practical calculations [HI, [l^. Due to the Pauli exclusion principle, fermionic systems have much 
higher energy scales than BECs. Consequently, a larger dipole moment is required to observe dipolar effects in dipolar 
Fermi gases than in dipolar BECs. Miyakawa et al. [l^ have shown that for heteronuclear molecules with a typical 
electric dipole moment of the order of 1 Debye, dipolar effects can be easily detected by studying the equilibrium 
properties of dipolar Fermi gases at zero temperature. On the other hand, Sogo et al. studied the dynamics of free 
expansion of dipolar Fermi gases [2^. They found that the aspect ratios of dipolar Fermi gases after free expansion 
are stretched along the dipole moment regardless of the trap geometry. This result contrasts to that of dipolar BECs, 
which strongly depends on the trap aspect ratio [Ij-Q- Some theoretical studies of dipolar Fermi gases have focused 
on the zero-temperature properties of these gases and there have been recent theoretical studies on finite-temperature 
gases [2i,[28H3i|. 

As described above, a huge dipole moment is required to observe dipolar effects in a Fermi system. However, 
despite many groups conducting experiments, no groups have succeeded in cooling polar molecules to the quantum- 
degenerate regime [6|. It is thus important to investigate the temperature range in which the dipolar interaction has 
appreciable effects. It is also important to quantify to what extent the dipolar interaction effect can be observed in 
the temperature regime of the current experiment. For this reason, we studied the equilibrium properties of dipolar 
Fermi gases above the Fermi temperature We have previously estimated the electric dipole moments and the 

trap frequencies required to observe deformations in momentum and real space in the high-temperature regime |3Q| . 
This result indicates that it is possible to observe deformation in real and momentum space in current experiments [6| . 

In the present paper, we study the equilibrium properties and stability of normal dipolar Bose gases above the Bose- 
Einstein transition temperature as well as dipolar Fermi gases. In particular, we focus on the different properties 
induced by the different statistics of dipolar Fermi and Bose gases because many experimental groups are progressing 
toward creating dipolar Bose molecules as well as dipolar Fermi molecules (e.g., ^^Rb^^K molecules [6| and ^^Rb'^^'^K). 
We show that the different properties of Bose and Fermi gases arise from their stabilities and the deformation of their 
momentum distributions. 

I. THEORY 

We consider Bose or Fermi heteronuclear molecules. The dipoles are assumed to be polarized along the z axis due 
to an external electric field. In the second quantized form, the Hamiltonian for this system is given by 



J dr^"^ (r) 



^(r) 



(1 + ^) 



g J dA"^ (r) (r) ^ (r) ^ (r) 



^ IJdrJ dr'^^ (r) (r^) y,, (r _ ^ (r^ ^ (r) , (1) 

where ^ (r) is the field operator and the hat indicates a second quantized operator. In the above equation, we defined 
77 = ±1, where we use the positive sign for the Bose case and the negative sign for the Fermi case. The first term in 
Eq. ([!]) describes the Hamiltonian of a single particle in a harmonic trap: 

V,.^,{r) = ^[o.l{x'+y')+coy], (2) 

where m is the particle mass. The second term describes the contact interaction due to 5-wave scattering, which 
occurs only for Bose gases, g is the coefficient of the interatomic interaction. The third term describes the two-body 
interaction Hamiltonian of the dipolar force: 

F,,(r) = ^(l-3^Vr^). (3) 

Here, d is the coefficient of the dipolar interaction, and d'^ = p^/47reo, where p is the magnitude of the electric dipole 
moment HH and eo is the electric permittivity of vacuum. 

To consider the phase-space distribution, we introduce the Wigner distribution function: 

W (p, r) = J dr'e-^P-^'/^(^^ (r - /2) ^ (r + r72)), (4) 



3 



The density distributions in real and momentum space in terms of (p, r) are given by n (p) = J dr W (p^r) and 
n (r) = {27rh)~ J dp VK(p,r), respectively. 

We consider the thermal equilibrium of a dipolar gas trapped in a harmonic potential. The system is not globally 
stable when a dipolar interaction is present because the interaction is partially attractive and causes collapse of the 
gas. However, a local metastable state exists at T = and at finite temperatures under certain conditions. As in 
Ref. [sO] , we look for this metastable state by introducing a variational Wigner distribution function. We assume the 
Maxwell-Boltzmann regime at relatively high temperatures, namely 



W (p, r) — exp ' 




Mo 



/ksT 



(5) 



where = 



= -\- y^^ = pI -\-py, and u = (cj^co'^)"^". Here, the positive parameters a and /3 respectively represent 
deformations of the density distributions in momentum and real space, and A describes isotropic compression in real 
space. These parameters are also introduced for T = in Ref. [l^- At finite temperatures, we introduce an additional 
variational parameter that characterizes isotropic compression in momentum space. The chemical potential is 
determined by the number constraint N = J n (r), which gives 
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From the normalization condition 



N 



dr 



dp 



(211%)' 



;W(p,r), 



(6) 



(7) 



the chemical potential /io is determined by Eq. (j6j). Under this ansatz, the density distribution in momentum space 
is given by 



2mkBT \ a 



(8) 



Thus, the aspect ratio of the momentum space distribution, which is the ratio of the root-mean-square mo mentum in 
the Pz direction to that in a direction in the Px-Py plane (we choose the Px direction in the following) \/ {vD / {pV)^ 
becomes Similarly, we obtain the density distribution in real space as 



n (r) = A/" 



2 \ 2 \ 3/2 



exp 



' 2kBT 



(9) 



leading to the aspect ratio of the real-space distribution (i.e., the ratio of the root-mean-square radius in the z direction 
to that in the x direction) being ^/ = We note that, from Eq. (j9j), A > 1 (A < 1) corresponds to 

compression (expansion) of the gas in real space, which is a consequence of the effective attraction (repulsion) of the 
dipolar interaction. In addition, comparing Eq. ([8]) with Eq. (|9j) reveals that plays the same role in real space as A 
plays in momentum space. 

To find a metastable state at a temperature T, we look for a local minimum of the Helmholtz free energy 



F = E-TS. 

Here, E is the total energy, and the entropy S is given by 

dp 



(10) 



-k> 



dp 



{W (p, r) In [W (p, r)] t[1±W (p, r)] In [1 ± (p, r)]} 



■W(p,r)ln[W(p,r)]. 



(11) 



In the first line of Eq. (fTTj) . the upper and lower signs are for Bose and Fermi gases, respectively. In the above 
expression for the entropy, we assumed W <C 1, which is consistent with the Maxwell-Boltzmann distribution. 

The total energy can be written as a sum of five contributions: the kinetic energy Ek^ the trapping potential energy 
Ey^ the interatomic mean-field energy Eg^ which arises only for Bose gases, the Hartree direct energy Eh-, and the 
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Fock exchange energy E^x- E = Ek + Ey + ^-^Eg + Eh + V^ex- The five contributions are given in terms of the 
Wigner distribution function as [30| 



(27r?i)'' J (27r?i)"' 

1 /■^j. /•^j., /• <ip /• dp' 

2 J ""J y (27r?i)^ J {2Trhf 



(14) 



Eh = - I dv I dv' I / -^Fdrf (r - r') W (p, r) (p', r') , (15) 



E.. = l/.R/.s/^/^y..(s)e<P-P')-V.^(p,K)^(p',R), 



(16) 



In Eq. ([Ej), we have introduced the center of mass coordinate R = (r + r') /2 and the relative coordinate s = r — r^ 
In deriving Eqs. ([15]) and (p!6]) . we used mean- field decoupling: 

(r) (rO ^ (r ) ^ (r)) 2^ (r) ^ (r))(^t (r^ ^ (r )) + (r) ^ (r ))(^'^ (r ) ^ (r)), (17) 

The difference of the sign of 77 in Eq. (p!6|) arises from the different statistics of the two systems. Consequently, the 
Fock exchange energy has different signs for Bose and Fermi systems. 

Inserting the variational function (|5]) into Eqs. (fT2]) - (fT6|) and ([TT]) . we express the free energy in terms of variational 
parameters: 



-3 + In N{XOuhf -3ln[kBT]^NkBT, (18) 

In the above equation, we have defined the scaling functions: 

/(a) = 2a + 1/^2, (19) 

gW) = 2po/P + l3V/3l (20) 

where Po = {ujp/ujz)^^'^ . The deformation function I {a) is defined as 

/ (a) ^ r de sin e ( ^ - 1) . (21) 



^Q;3 sin^ G + cos^ 
The property of / (a) is described in Appendix [A] In addition, we defined 

in Eq. (p!8|) . The first and second terms in the first line of Eq. (fTSj) are the kinetic and potential energies, respectively. 
The kinetic energy is inversely proportional to the square of d. Thus, an increase in d indicates a reduction in the 
effective temperature T* = T jd^ associated with the kinetic energy. The third term is the interatomic mean- field 
energy due to 5- wave scattering. The interatomic mean- field energy contributes only to the compression/expansion 
in the real-space distribution. Although the mean-field energy is isotropic, it is closely related to the stability of 
the Bose system, as discussed below. The fourth and fifth terms are respectively the Hartree direct energy and the 
Fock exchange energy due to the dipolar interaction. Note that the deformation parameters a and /3 contribute 
independently to the internal energy, E. The momentum space deformation parameter a appears only in the kinetic 
and Fock exchange energies, whereas the real space deformation parameter /3 appears only in the potential and Hartree 
direct energy. The last three terms in Eq. (p!8|) express the entropy, which is independent of a and /3. 
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II. EQUILIBRIUM SOLUTIONS 

We now discuss equilibrium solutions of dipolar quantum gases. The equilibrium solution can be found by mini- 
mizing the free energy with respect to the four variational parameters, a, /3, A, and 9. To be a local minimum of the 
free energy representing a met ast able state, the solution must be a convex downward point in the four-dimensional 
space of the variational parameters. This requires the following conditions: 

The above conditions arise from dF/da = and dF/df3 = exploiting the properties of f (a) and g{/3). The 
difference between Fermi and Bose systems arise from the different signs of the Fock exchange energy. 

From the conditions (|23|) . we find that the momentum space distribution is always elongated in the pz direction in 
a Fermi gas, whereas it is stretched in the Px-Py plane in a Bose gas. In contrast, the real-space distribution tends to 
be stretched in the z direction and the aspect ratio of the real-space distribution is larger than that for non-dipolar 
gases in both Fermi and Bose systems. 

These different properties can be easily explained. If there is no Fock exchange energy, the conditional equations, 
a = 1 and ^ = 1, apply for both Fermi and Bose systems. The Fock exchange energy leads to deformation and 
compression of the momentum distribution. Because the sign of the Fock exchange energy differs for Bose and Fermi 
systems, deformation in momentum space occurs in opposite directions in the two systems. 

To find the equilibrium solution, we numerically minimize the free energy by varying three system parameters 
(temperature T, trap aspect ratio /3o, and the magnitude of the electric dipole moment p or the ratio between the 
dipolar and interatomic interactions g/(P) for m = 100 a.u.m., = 27r x 10^ Hz, and TV = 10^. We normalize the 
temperature using the ideal gas Fermi temperature = {6N)^^^ Tiuj/kB — 188 nK for a Fermi system and the ideal 
gas BEC temperature = [TiujN^^'^^ / (1.202-*^/^/^^) — 97.2 nK for a Bose system, which are defined for an isotropic 
trap. 

A. The interaction energy dependencies of the variational parameters 

First, we analyze the variational parameters as functions of the interaction energy. Figure [1] shows plots of the 
variational parameters of the dipolar Fermi gas as functions of the dipole moment for different trap aspect ratios 
when the temperature is fixed at T= 2Tp. For large trap aspect ratios, we plotted the optimized values for p < Pc^ 
where Pc is the critical electric dipole moment above which the variational free energy (fT8|) has no local minimum 
and the system starts to collapse. For example, the critical dipole moment is Pc = 4.36 Debye for /3o = 1.0 and 
Pc = 3.30 Debye for /3o = 2.0. On the other hand, systems with /3o = 0.5 are stable for p < 6.0 Debye. These 
results reveal that the system is unstable in highly elongated cigar-shaped traps and that the system tends to become 
unstable with increasing dipolar interaction. Figure [T]^a) shows that the momentum distribution is elongated in the 
dipole direction regardless of the trap aspect ratio. It is important to note that no deformation of the anisotropic 
momentum distribution is induced without the dipolar interaction. Next, Fig.djb) shows that with increasing dipolar 
interaction, the real-space distribution becomes more elongated in the z direction relative to the real-space distribution 
without the dipolar interaction. In addition, an isotropic trap gives the largest shift in /3. Next, Fig.[T]^c) reveals that 
A changes from A < 1 to A > 1 with increasing trap aspect ratio. The reason for this change is that as the trap shape 
changes from being oblate to being cigar-shaped, the dipolar interaction changes from being repulsive to attractive 
due to anisotropy of the dipolar interaction. Finally, from Fig. [TJd), we find that the effective temperature decreases 
as the dipolar interaction increases regardless of the trap aspect ratio because only the kinetic energy and entropy 
depend on 6 (Eq. (ITSl) ). 

Figure [2] shows plots of the scaling parameters of the dipolar Bose gas as functions of the ratio of the interatomic 
to dipolar interactions g/(P for different trap aspect ratios with the temperature fixed to T= 2T^ . Systems with 
— 0.5, /3 = 1.0, and /3 = 2.0 are stable in the region —2 < {g/d?) < 6. Figure [2j a) shows that the momentum 
distribution is stretched perpendicular to the dipolar direction contrary to the case of the dipolar Fermi gases. In 
contrast. Fig. [2fb) reveals that with increasing effective dipolar interaction, the real-space distribution becomes more 
elongated in the z direction relative to the case without the dipolar interaction and that j3 is shifted the most in an 
isotropic trap, just as for dipolar Fermi gases. Next, Fig. [2fc) shows that with decreasing interaction ratio, A shifts 
to A > 1. This is because the large positive interaction ratio indicates that the interatomic interaction is large and 
repulsive, which expands the distribution (A < 1). On the other hand, the large negative interaction ratio means that 
the interatomic interaction is large and attractive, which compresses the distribution (A > 1). In addition, the density 
increases with decreasing interaction ratio, which enhances the effective dipolar interaction. Finally, from Fig. [2fd), 
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we find that the effective temperature decreases as the dipolar interaction increases regardless of the trap aspect ratio 
because only the kinetic energy and entropy depend on (Eq. ([lEj))- 

Comparing Figs. [T] and [21 we find that the different statistics of fermions and bosons manifests itself only as 
deformation in momentum space and compression in real space. First, dipolar Fermi and Bose gases have opposite 
properties in momentum space because they have opposite signs for the Fock exchange energy. Next, because only 
Bose systems have an interatomic interaction, which contributes to real-space compression, dipolar Bose gases depend 
on the sign of the ratio of the interatomic interaction to the dipolar interaction in addition to their strengths. The 
interatomic interaction has a significant effect on the stability, which we show later. 



B. The temperature dependence of the variational parameters 

1. The temperature dependence of the variational parameters for several trap aspect ratios 

Next, we analyze the variational parameters as a function of temperature. Figure [3] shows a plot of the variational 
parameters of the dipolar Fermi gas as functions of temperature for j3o = 0.4 (solid line), j3o = 1.0 (dashed line), 
and /3o = 2.0 (dash-dotted hue), where the electric dipole moment is fixed at p = 1 Debye. Figure [3f a) indicates 
that the momentum distribution becomes more elongated in the dipole direction as the temperature decreases. The 
deformation in the momentum distribution arises from competition between the kinetic energy, which favors an 
isotropic momentum distribution, and the exchange energy, which favors an anisotropic momentum distribution. 
With decreasing temperature, the kinetic energy decreases and the interaction energy increases, increasing the aspect 
ratio in the momentum distribution. Figure O^b) reveals that a lower temperature increases the deviation of the 
aspect ratio in real-space from that of the non-interacting case. It also shows that the isotropic trap has the largest 
shift of p. Figure [3l^c) shows that the magnitude of the deviation from 1 in A increases at lower temperatures. In 
addition, we have A > 1 for /3o > 1 (cigar-shaped trap) and A < 1 for /3o < 1 (oblate trap). The Hartree direct energy 
becomes attractive (repulsive) in the former (latter) case due to deformation of the spatial density distribution, which 
is mainly determined by the trap geometry. This is leads to compression (expansion) of the gas cloud for /3o > 1 
(/3o < !)• Finally, Fig. [3fd) reveals that the dipolar Fermi gas is compressed in momentum space irrespective of the 
trap aspect ratio. The interaction effect increases with decreasing temperature in all cases, whereas the deformations 
in real and momentum spaces are negligibly small for the case T > 2T^, p = 1 Debye, and a; = 27r x 10^ Hz. 
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FIG. 2: Variational parameters of the dipolar Bose gas as functions of the ratio of the interatomic interaction to the dipolar 
interaction g/dP at the temperature T= 2.0Tc for a trap asymmetry /3o = 0.5 (solid line), /3o = 1.0 (dashed line), and /3o = 2.0 
(dash-dotted line). 



Figure [4] shows a plot of the variational parameters of the dipolar Bose gas as functions of temperature for /3o = 0.4 
(solid line), /3o = 1.0 (dashed line), and /3o = 2.0 (dash-dotted line) and for when the interaction ratio is fixed at 
{g/dP) = 1.0. Figure ID^a) indicates that the momentum distribution becomes more stretched perpendicular to the 
dipole direction as the temperature decreases, which is the contrary to the Fermi case. Figure [H^b) reveals that the 
deformation in real space is the same as that for a dipolar Fermi gas. Figure [4](c) shows that the magnitude of 
the deviation in A increases at lower temperatures because the effective dipolar interaction increases with decreasing 
temperature. As the dipolar interaction depends on the trap aspect ratio, we have A > 1 for /3o > 1 and A < 1 for 
/3o < 1. Finally, Fig. HJ^d) reveals that a dipolar Fermi gas is compressed in momentum space irrespective of the trap 
aspect ratio. 



2. The temperature dependence of the variational parameters for different interaction strengths 

Figure [5] shows the variational parameters of the dipolar Fermi gas as functions of temperature for p = 1.0 Debye 
(solid line), p = 3.0 Debye (dashed line), and p = 5.0 Debye (dash-dotted line) for an isotropic trap /3o = 1.0. In all 
cases, the deviation of the variational parameters from those of the non-interacting case becomes more pronounced at 
lower temperatures and at higher electric dipole moments. For large electric dipole moments, we plotted the optimized 
values for T > Tc. The critical temperature for this collapse for p = 5.0 Debye is Tc = 2.25Tp. On the other hand, 
systems with p = 1.0 Debye and p = 3.0 Debye are stable in the temperature regime T > 1.5Tp. These results reveal 
that the system is unstable even in the high-temperature regime for large electric dipole moments. Figures [5ja) and 
(b) show that the aspect ratios in both momentum and real spaces are sufficiently large to observe deformation effects 
in experiments for polar molecules with large electric dipole moments. 

Figure [6] shows the variational parameters of the dipolar Bose gas as a function of temperature for g/dP = —1.0 
(solid line), g/cP = 0.0 (dashed line), and g/cP = 1.0 (dash-dotted line) for an isotropic trap /3o = 1.0. In all cases, 
the deviation of the variational parameters from those for the non-interacting case becomes more pronounced at 
lower temperatures and at smaller ratios of the interactions. This is because a small interaction ratio indicates a 
small repulsive or attractive interatomic interaction, which shortens the atomic distance and enhances the dipolar 
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FIG. 3: Variational parameters of the dipolar Fermi gas as functions of the temperature T for a fixed dipole moment p = 1 Debye 
and for trap asymmetries of f3o — 0.4 (solid line), /3o = 1.0 (dashed line), and /3o = 2.0 (dash-dotted line). 



interaction. We plotted the optimized values for T > Tc for a small interaction ratio. The critical temperature for this 
collapse is = IMT^ for g/d'^ = -1.0. On the other hand, g/d^ = 0.0 and g/d^ = 1.0 are stable in the temperature 
regime T > 1.5T^. In addition, from Fig. [6l we see that the magnitude of the deviation from 1 in A depends on the 
sign of the interaction ratio at low temperatures. 



III. STABILITY 
A. Fermi gas 

In this section, we discuss the stability of dipolar gases. We first consider a Fermi gas. Figure [7] shows the stability 
diagram in which the critical values of the electric dipole moment Pc are plotted as functions of the trap aspect ratio 
for different temperatures. This result indicates that the critical dipole moment pc increases drastically as the 
trap becomes more oblate. This is consistent with the result for a zero-temperature Fermi gas [l^- The unstable 
region expands with decreasing temperature. Based on Fig. [5l we assume that the critical temperature for collapse 
instability increases as the electric dipole moment increases. This behavior is confirmed by Fig.[8l^a), which shows the 
critical temperature Tc for f3o = 0.5, /3o = 1.0, and /3o = 2.0. Figure [Hl^b) shows the critical temperature as a function 
of the trap aspect ratio for p = 3.0 and p = 5.0 Debye. This result indicates that the critical temperature Tc increases 
as /3o increases. When p= 1.0 Debye, the system is always stable within the parameter range shown in the figure 
{Po < 40). Figure [Hfb) suggests that we should be careful not to enter the unstable region when the temperature is 
reduced. Even if the system is initially in a stable region at a high temperature, the system may (depending on the 
trap aspect ratio) become unstable on cooling before reaching the quantum degenerate regime. In addition, the local 
minimum always disappears when A exceeds a critical value. This suggests that the instability first occurs in real 
space and a dipolar Fermi gas collapses when A reaches the critical value. 
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FIG. 4: Variational parameters of the dipolar Bose gas as functions of the temperature T for a fixed ratio between the dipolar 
and interatomic interactions g/dP = 1.0 for trap asymmetries of /3o = 0.4 (solid line), /3o = 1.0 (dashed line), and /3o = 2.0 
(dash-dotted line). 



B. Bose gas 



Next, we consider Bose gases. Figure [9] shows the stability diagram in which the critical values of the ratio of 
the interatomic interaction to the dipolar interaction (^g/dP^ ^ are plotted as functions of the trap aspect ratio /3o for 
different temperatures. As the interaction ratio approaches 0, the dipolar interaction is enhanced, so that the stability 
becomes more dependent on the trap geometry. Increasing the attractive interatomic interaction causes the system 
to become unstable. Meanwhile, the effective repulsive dipolar interaction in the oblate trap stabilizes the system. 
This means that the dipolar interaction keeps the system stable against the attractive interatomic interaction at low 
trap aspect ratios (/3o < 1). This result is consistent with Fig. [2l In Fig. [TOT a). we plot the critical temperatures 
as functions of the ratio g/d? for j3o = 1.0 and /3o = 2.0. The dependence of the trap aspect ratio on the stability 
decreases with decreasing interaction ratio because the instability is caused by the attractive interatomic interaction, 
and not by the dipolar interaction. When /3o = 0.5, the system is always stable within the parameter range shown 
in the figure (—8 < g/dP < 4). Figure [TOTb) shows the critical temperature as functions of the trap aspect ratio for 
different ratios of the interactions. This result indicates that the critical temperature increases with increasing trap 
aspect ratio. In Fig. [TOTb). the critical temperatures only shift down with increasing ratio of the interaction because 
expansion or compression of the distribution arising from the interatomic interaction is independent of the trap aspect 
ratio. In addition, we see that the local minimum always disappears when A exceeds a critical value, as in the case of 
Fermi gases. This suggests that the instability first occurs in real space for both Fermi and Bose gases. 



IV. RELATION TO EXPERIMENTS 



We discuss the stability properties in connection with experiments. First, we consider the Fermi case. The JILA 
experiment Q revealed that the mass of the polar molecule is about m= 128 a.m.u. and that the electric dipole 
moment is about p= 0.566 Debye. Taking these values along with uo = 27r x 10^ and = 10^, we find that it is 
always possible to find a stable dipolar Fermi gas for any temperature T > 1.5T^ in the trap aspect ratio regime 
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c) 



d) 





FIG. 5: Variational parameters of the dipolar Fermi gas as functions of the temperature for a fixed trap asymmetry of /3o = 1.0 
for a dipole moment p — 1.0 Debye (solid line), p = 3.0 Debye (dashed line), and p = 5.0 Debye (dash-dotted line). 



/3o < 40. Since, according to the results in Ref. [17], a gas of polar molecules at zero temperature with the same 
system parameters becomes unstable, the critical temperature should lie in < T < 1.5Tp. It might be expected 
that the deformation effect of the distribution can be observed using a highly elongated cigar-shaped trap potential. 
However, the dipolar Fermi gas becomes unstable at extremely large trap aspect ratios. For uj = 27r x 10^, N = 10^, 
p = 0.566 Debye, and m = 128 a.m.u., the critical trap aspect ratios for T = 1.5T^ and T = 2.0T^ are /^g = 62.0 and 
/3q = 72.5, respectively. These effects suggest that deformation effects cannot be observed by just adjusting the trap 
aspect ratio. 

We now discuss how the deformation effects of momentum-space and real-space distributions can be observed in 
the polar molecules used in the JILA experiment [6]. We assume that p, and m have the same values as those 
given in the previous paragraph. Increasing u causes dipolar interaction effects to become more pronounced because 
the density of the gas trapped in a tight potential increases and the interaction becomes stronger. For = 57r x 10^ 
Hz, /3o= 1.0, and T = 1.5Tp, we have aspect ratios \/ {pI)/{Px) — <^~^^^ — 1.130 in the momentum distribution and 
y^{z'^)/{x^) = 1.127 in the spatial distribution. For uj = 67t x 10^ Hz, /3o= 1.0 and T = 2.0T^, we have 

^-3/2 ^ 1.106 and 1.103. Thus, the deformation effects arising from the anisotropic nature of the dipolar 

interaction are detectable using the present experimental conditions. 

Next, we consider the Bose case. We assume ^^K^^Rb molecule 0, so we cite the physical quantities from the JILA 
experiment [6^, which uses isomeric molecules: m= 128 a.m.u. and p= 0.566 Debye. Taking these values along with 
a scattering length of a = 10 nm, a coefficient of the dipolar interaction of (P ^ 3.20 x 10~^^, and an interatomic 
interaction of ^ 6.58 x 10~^^, the ratio of those interactions is g/(P ~ 0.21. Let us now consider the possibility of 
experimentally observing the deformation caused by the dipolar interaction. First, we consider a highly elongated trap 
potential with u = 27r x 10^ and N = 10^. At the temperature T = 2.0T^, the system is stable for (3o = 100.0 but the 
deformation in momentum and real spaces scaled by non-dipolar gases are only ^ Q (/3//3o)^^^ - 1.000. 

On the other hand, at a temperature of T = 1.5T^, the aspect ratios are ot'^I^ :^ 0.959 and (/3//3o)^^^ — 1.000 for 
this trap aspect ratio. These results indicate that using a highly elongated cigar-shaped trap potential has limitations 
and is ineffective, especially in real space. 

We next consider large trap potentials using the same values for TV, gj^^ and m as those given in the previous 
paragraph. For a; = IItt x 10^ Hz, /3o= 1.0, and T = 1.5T^^, we have aspect ratios ot'^I^ ~ 0.912 in the momentum 
distribution and /^^/^ 1.091 in the spatial distribution. For a; = 307r x 10^ Hz, /3o= 1.0 and T = 2.0T^°, we have 
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FIG. 6: Variational parameters of the dipolar Bose gas as functions of the temperature for a fixed trap asymmetry /3o ■ 
for a ratio between the dipolar and interatomic interactions g/(P = —1.0 (solid line), g/dP — 0.0 (dashed line), and g/dP 
(dash-dotted line). 



1.0 
1.0 




FIG. 7: Critical value of the dipole moment pc of the dipolar Fermi gas as a function of the trap aspect ratio for T /Tp = 1.5 
(open circles), T/T^ = 2.0 (filled circles), and T/T^ = 3.0 (crosses). 



a ~ 0.930 and ^ 1.072. Thus, the deformation effects can be observed by using a high trap frequency rather 
than by using a high trap aspect ratio. 



V. CONCLUSION 



We have used a variational method to study the equilibrium properties of a dipolar gas at finite temperatures. For 
a Fermi gas at zero temperature, the anisotropic nature of the dipolar interaction leads to deformations in momentum 
and real space, and the partial attraction of the interaction causes the gas to collapse. In addition, the dipolar Fermi 
gas becomes compressed in momentum space as the electric dipole moment increases. We found that the deformations 
in both momentum and real space can be observed in the high-temperature regime with a large electric dipole moment 
and a high trap frequency. In addition, the stable region expands at finite temperatures. These results will be useful 
when cooling polar molecules in experiments. 

By studying a Bose gas, we find that the different statistics of Fermi and Bose gases gives rise to different signs for 
the Fock exchange energy as well as the emergence of the interatomic interaction in dipolar Bose gases. The different 
signs of the Fock exchange energy causes opposite properties: the momentum distribution is elongated in the dipole 
direction in a dipolar Fermi gas, whereas the momentum distribution is stretched perpendicular to the dipole direction 
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P Po 

FIG. 8: Critical temperature of the dipolar Fermi gas. (a) Critical temperature as a function of the dipole moment for trap 
aspect ratios of /3o = 0.5 (open circles), /3o = 1.0 (filled circles), and /3o = 2.0 (crosses), (b) Critical temperature as a function 
of the trap aspect ratio for dipole moments of p = 3.0 (filled circles) and p = 5.0 (crosses). This small dip between /3o = 20 and 
28 for p = 5.0 is due to numerical uncertainty. 




FIG. 9: Critical value of the ratio between the dipolar interaction and the mean- field energy [g/dP^ ^ of the dipolar Bose gas 
as a function of the trap aspect ratio for T/T^ — 1.5 (open circles), T/T^ — 2.0 (filled circles), and T/T^ — 3.0 (crosses). 



in a dipolar Bose gas. 

In addition, we found that the local minimum of the free energy always disappears when A exceeds a critical value 
for both Fermi and Bose gases. This implies that in both cases the instability first occurs in real space. However, 
the different statistics mainly arise from the sign of the Fock exchange energy, which indicates that the difference 
mainly arises in momentum space. Consequently, there are few differences in the critical interaction and the critical 
temperature as a function of the trap aspect ratio for Fermi and Bose gases. On the other hand, there is a clear 
difference in the critical temperature as a function of the interaction for the two systems. The dipolar interaction 




FIG. 10: Critical temperature of the dipolar Bose gas. (a) Critical temperature as a function of the ratio between the dipolar 
interaction and the mean-field energy for trap aspect ratios /3o = 1.0 (filled circles) and /3o = 2.0 (crosses), (b) Critical 
temperature as a function of the trap aspect ratio for the ratio between the dipolar interaction and the mean-field energy 
g/d? = —1.0 (open circle), g/dP = 0.0 (filled circles), and g/d^ = 1.0 (crosses). 
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changes from being attractive to being repulsive on decreasing the trap aspect ratio. On the other hand, the repulsive 
interatomic interaction keeps the system stable, and the attractive interatomic interaction causes the system to become 
unstable. For this reason, decreasing the ratio of the interatomic interaction to the dipolar interaction causes the 
instability to become independent of the trap aspect ratio for a Bose gas. 

Our variational formalism can be extended to study the dynamics of dipolar gases at finite temperatures, such as 
the expansion dynamics and collective oscillations. In particular, expansion dynamics is important since it can provide 
direct information about deformation of the momentum distribution. This type of dynamics is also useful in studying 
the different statistics, which is one of the main results of this study. In addition, it may be interesting to study the 
expansion dynamics of a dipolar Bose-Fermi mixture. We are also interested in the equilibrium and dynamics below 
the Bose-Einstein transition temperature for dipolar Bose gases. We note that the relation between the real-space 
distribution and the momentum-space distribution of BECs is described by the Fourier transform of the macroscopic 
wave function, but it is not the case for a thermal cloud gas. Therefore, fascinating properties are expected below the 
Bose-Einstein transition temperature. 

In the future, we intend to investigate the temperature dependence of the expansion dynamics and collective mode 
in the presence of the dipolar interaction. In addition, we will study the equilibrium and dynamics of a dipolar 
Bose-Fermi mixture and dipolar Bose gases below the Bose-Einstein transition temperature. We hope this study will 
stimulate further experiments on dipolar gases at finite temperatures. 
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Appendix A: Deformation function 



In this appendix, we consider the deformation function (|2T]) . which determines the properties of the Fock exchange 
energy and the Hartree direct energy. Rewriting the deformation function, one can derive the analytical expression |32| : 



J (a) 



aB sm B I — ^ ^ — 

V sin B + cos^ B 



1 




(a<l) 
(a=l) 
{a > 1) 



(Al) 



This formula indicates that / (a) is a monotonically decreasing function of a with — 2 < / (a) < 4; it is positive for 
a < 1, passes through zero at a = 1, and is negative for a > 1 Eli- Figure [TT] shows the deformation function 
and its derivative. 




dl(a)/da 




FIG. 11: Deformation function I (a) and the derivative of I (a) as a function of a. 
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